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The three-dimensional topological insulator (originally called "topological insulators") is the first 
example in nature of a topologically ordered electronic phase existing in three dimensions that can- 
not be reduced to multiple copies of quantum-Hall-like states. Their topological order can be 
realized at room temperatures without magnetic fields and they can be turned into 
magnets and exotic superconductors leading to world-wide interest and activity in 
topological insulators. These experimentally demonstrated unique properties of 3D topological 
insulators have turned the topic into a field of intense experimental activity. It is also only the 
third topologically ordered phase in weakly interacting systems to be discovered in nature, with 
the first two belonging to the quantum Hall-like topological insulator class consisting of the 2D 
integer quantum Hall state (IQH) and the 2D quantum spin Hall state (QSH). The 2D quantum 
spin Hall state (QSH) can be thought of as two copies of IQH put together leading to a time- 
reversal invariant version of IQH state. All of the 2D topological insulator examples (IQH, QSH) 
including the fractional one (FQH) involving Coulomb interaction are understood in the standard 
picture of quantized electron orbits in a spin-independent or spin-dependent magnetic field, the 3D 
topological insulator defies such description and is a novel type of topological order 
which cannot be reduced to multiple copies of quantum-Hall-like states. In fact, the 3D 
topological insulator exists not only in zero magnetic field, and differs from the 2D variety in three 
very important aspects: 1) they possess topologically protected 2D metallic surfaces (a new type of 
2DEG) rather than the ID edges, 2) they can work at room temperature (300K and beyond, large- 
gap topological insulators) rather than cryogenic (mK) temperatures required for the QSH effects 
and 3) they occur in standard bulk semiconductors rather than at buried interfaces of ultraclean 
semiconductor heterostructures thus tolerate stronger disorder than the IQH-like states. One of 
the major challenges in going from quantum Hall-like 2D states to 3D topological insulators is to 
develop new experimental approaches/methods to precisely probe this novel form of 
topological-order since the standard tools and settings that work for IQH-state also work for QSH 
states. The method to probe 2D topological-order is exclusively with charge transport (pioneered 
by Von Klitzing in the 1980s), which either measures quantized transverse conductance plateaus in 
IQH systems or longitudinal conductance in quantum spin Hall (QSH) systems. In a 3D topological 
insulator, the boundary itself supports a two dimensional electron gas (2DEG) and transport is not 
(Z2) topologically quantized hence cannot directly probe the topological invariants Vo or the topo- 
logical quantum numbers analogous to the Chern numbers of the IQH systems. This is unrelated 
to the fact that the present materials have some extrinsic or residual/impurity conductivity in their 
naturally grown bulk. In this paper, we review the birth of momentum- and spin-resolved spec- 
troscopy as a new experimental approach and as a directly boundary sensitive method to study and 
prove topological-order in three-dimensions via the direct measurements of the topological invariants 
fo that are associated with the Z2 topology of the spin-orbit band structure and opposite parity 
band inversions, which led to the experimental discovery of the first 3D topological insulator in Bi- 
Sb semiconductors (KITP Proceeding 2007 http://online.itp.ucsb.edu/online/motterials07/hasan/ 1 
(2007), Nature 452, 970 (2008), Submitted 2007) which further led to the discovery of the Bi2Se3 
class - the most widely researched topological insulator to this date. We discuss the fundamental 
properties of the novel topologically spin-momentum locked half Dirac metal on the surfaces of the 
3D topological insulators and how they emerge from topological phase transitions due to increasing 
spin-orbit coupling in the bulk. These electronic and spin properties of topological surface states 
discovered via the methods reviewed here are already guiding the interpretation of surface transport 
measurements as they are beginning to be possible further advancing the field potentially towards 
device applications. These methods and their derivatives are also being applied by many others 
world-wide for further investigations of topological-order and for discovering new topological insula- 
tor states as well as exotic topological quantum phenomena (the list is too long to review here) . We 
also review how spectroscopic methods are leading to the identification of spin-orbit superconduc- 
tors that may work as Majorana platforms and can be used to identify topological superconductors 
- yet another class of new state of matter. 
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THE BIRTH OF MOMENTUM-RESOLVED 
SPECTROSCOPY AS A DIRECT 
EXPERIMENTAL PROBE OF 
TOPOLOGICAL-ORDER 

Ordered phases of matter such as a superfluid or a 
ferromagnet are usually associated with the breaking of 
a symmetry and are characterized by a local order pa- 
rameter [1]. The typical experimental probes of these 
systems are sensitive to order parameters. In the 1980s, 
two new phases of matter were realized by subjecting 2D 
electron gases at buried interfaces of semiconductor het- 
erostructures to large magnetic fields. These new phases 
of matter, the 2D integer and 2D fractional quantum Hall 
states, exhibited a new and rare type of order that is 
derived from an organized collective quantum entangled 
motion of electrons These so-called "2D topolog- 

ically ordered insulators" do not exhibit any symmetry 
breaking and are characterized by a topological number 
[il as opposed to a local order parameter. The most 
striking manifestation of this 2D topological order is the 
existence of one-way propagating ID metallic states con- 
fined to their edges, which lead to remarkable quantized 
charge transport phenomena. To date the experimen- 
tal probe of their topological quantum numbers is based 
on charge transport, where measurements of the quanti- 
zation of transverse magneto-conductivity axy = ne'^/h 
(where e is the electric charge and h is Planck's constant) 
reveals the value of the topological number n that char- 
acterizes these quantum Hall states [7]. 

Recently, a third type of 2D topological insulator, the 
spin quantum Hall insulator, was theoretically predicted 
[H, [il and then experimentally discovered [13]. This class 
of quantum Hall-like topological phases can exist in spin- 
orbit materials without external magnetic fields, and can 
be described as an ordinary quantum Hall state in a spin- 
dependent magnetic field. Their topological order gives 
rise to counter-propagating ID edge states that carry op- 
posite spin polarization, often described as a superposi- 
tion of a spin up and spin down quantum Hall edge state 
(Figure 1). Like conventional quantum Hall systems, the 
2D spin quantum Hall insulator (QSH) is realized at a 
buried solid interface. The first, and to date only, real- 
ization of this phase was made in (Hg,Cd)Te quantum 
wells using charge transport by measuring a longitudinal 
conductance of about 2e^ /h at mK temperatures [isj. 
The quantum spin Hall state (QSH) can be thought of 
as two copies of integer quantum Hall states (IQH) and 
protected by a Z2 invariant. 

It was also realized that a fundamentally new type of 
genuinely three-dimensional topological order might be 
realized in bulk crystals without need for an external 
magnetic field p^Ql-[l2|. Such a 3D topological insulator 
cannot be reduced to multiple copies of the IQH and such 
phases would be only the fourth type of topologically 
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FIG. 1: 2D and 3D topological insulators in nature. The 2D topological insulator class consists of the 2D quantum Hall 
states (IQH) and 2D quantum spin Hall state (QSH). The latter is constructed from two copies of the former. On the other 
hand, a 3D quantum Hall state is forbidden in nature, so the 3D topological insulator represents a new type of topologically 
ordered phase. The protected surface states form a novel type of topological metal (half Dirac metal) where electron's spin is 
locked to its momentum but exhibit no spin quantum Hall effect. 



ordered phase to be discovered in nature, and the first 
type to fall outside the quantum Hall-like 2D topological 
states (IQH, FQH, QSH). Instead of having quantum- 
Hall type ID edge state, these so-called 3D topological in- 
sulators would possess unconventional metallic 2D topo- 
logical surface states called spin-textured helical metals, 
a type of 2D electron gas long thought to be impossible to 
realize. However, it was recognized that 3D topological 
insulators would NOT necessarily exhibit a topologically 
(Z2) quantized charge transport by themselves as carried 
out in a conventional transport settings of all quantum- 
Hall-like measurements. Therefore, their 3D topological 
quantum numbers (Z2), the analogues of n (Chern num- 
bers), could not be measured via the charge transport 
based methods even if a complete isolation of surface 
charge transport becomes routinely possible. Owing to 
the 2D nature of the two surface conduction channels that 
contribute together in a 3D topological insulator case, it 
was theoretically recognized that it would not be possi- 
ble to measure the topological invariants due to the lack 
of a quantized transport response of the 2D surface that 



measures the Z2 topological invariants [12|. 

Here we review the development of spin- and angle- 
resolved photoemission spectroscopy (spin-ARPES) as 
the new approach/ method to probe 3D topological or- 
der [13, [ill, which today constitutes the experimental 
standard for identifying topological order in bulk solids 
also used by many others world-wide. We will review the 
procedures for i) separating intrinsic bulk bands from 
surface electronic structures using incident energy mod- 
ulated ARPES, ii) mapping the surface electronic struc- 
ture across the Kramers momenta to establish the topo- 
logically non-trivial nature of the surface states, iii) us- 
ing spin-ARPES to map the spin texture of the sur- 
face states to reveal topological quantum numbers and 
Berry's phases and iv) measuring the topological par- 
ent compounds to establish the microscopic origins of 3D 
topological order. These will be discussed in the context 
of Bii_a^Sba^, which was the first 3D topological insulator 
to be experimentally discovered in nature and a textbook 
example of how this method is applied. The confluence of 
three factors, having a detailed spectroscopic procedure 
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to measure 3D topological order, their discovery in fairly 
simple bulk semiconductors and being able to work at 
room temperatures, has led to worldwide efforts to study 
3D topological physics and led to over 100 compounds 
being identified as 3D topological insulators to date. 

SEPARATION OF INSULATING BULK FROM 

METALLIC SURFACE STATES USING 
INCIDENT PHOTON ENERGY MODULATED 
ARPES 

Three-dimensional topological order is predicted to oc- 
cur in semiconductors with an inverted band gap, there- 
fore 3D topological insulators are often searched for in 
systems where a band gap inversion is known to take 
place as a function of some control parameter. The ex- 
perimental signature of being in the vicinity of a bulk 
band inversion is that the bulk band dispersion should 
be described by the massive Dirac equation rather than 
Schroedinger equation, since the system must be de- 
scribed by the massless Dirac equation exactly at the 
bulk band inversion point. 

The early theoretical treatments ^ol, Hil focused on 
the strongly spin-orbit coupled, band- inverted Bii_a^Sba^ 
series as a possible realization of 3D topological or- 
der for the following reason. Bismuth is a semimetal 
with strong spin-orbit interactions. Its band structure 
is believed to feature an indirect negative gap between 
the valence band maximum at the T point of the bulk 
Brillouin zone (BZ) and the conduction band minima 
at three equivalent L points [l^, [HI (here we gener- 
ally refer to these as a single point, L). The valence 
and conduction bands at L are derived from antisym- 
metric (La) and symmetric (L^) p-type orbitals, respec- 
tively, and the effective low-energy Hamiltonian at this 
point is described by the (3+l)-dimensional relativistic 
Dirac equation [26- 28]. The resulting dispersion relation, 

E{k) = :^\J (y -k) -f ^ v ■ \s highly linear owing 
to the combination of an unusually large band velocity 
V and a small gap A (such that |A/|iT|| ^ 5 x 10~^A~^) 
and has been used to explain various peculiar properties 
of bismuth [26l-[28| . Substituting bismuth with antimony 
is believed to change the critical energies of the band 
structure as follows (see Figl2]). At an Sb concentration 
of X ~ 4%, the gap A between and closes and a 
massless three-dimensional (3D) Dirac point is realized. 
As X is further increased this gap re-opens with inverted 
symmetry ordering, which leads to a change in sign of A 
at each of the three equivalent L points in the BZ. For 
concentrations greater than x ^ 7% there is no overlap 
between the valence band at T and the conduction band 
at L, and the material becomes an inverted-band insula- 
tor. Once the band at T drops below the valence band 
at L, at X ^ 8%, the system evolves into a direct-gap 
insulator whose low energy physics is dominated by the 
spin-orbit coupled Dirac particles at L [10, 24]. 

High-momentum-resolution angle-resolved photoemis- 



sion spectroscopy performed with varying incident pho- 
ton energy (IPEM-ARPES) allows for measurement 
of electronic band dispersion along various momentum 
space (/c-space) trajectories in the 3D bulk BZ. ARPES 
spectra taken along two orthogonal cuts through the L 
point of the bulk BZ of Bio.gSbo.i are shown in Figs la 
and c. A A-shaped dispersion whose tip lies less than 50 
meV below the Fermi energy {Ep) can be seen along both 
directions. Additional features originating from surface 
states that do not disperse with incident photon energy 
are also seen. Owing to the finite intensity between the 
bulk and surface states, the exact binding energy (Eb) 
where the tip of the A-shaped band dispersion lies is un- 
resolved. The linearity of the bulk A-shaped bands is 
observed by locating the peak positions at higher Eb 
in the momentum distribution curves (MDCs), and the 
energy at which these peaks merge is obtained by ex- 
trapolating linear fits to the MDCs. Therefore 50 meV 
represents a lower bound on the energy gap A between 
La and Lg. The magnitude of the extracted band veloc- 
ities along the kx and ky directions are 7.9±0.5 x 10^ 
ms~^ and 10.0 ± 0.5 x 10^ ms~^, respectively, which are 
similar to the tight binding values 7.6 x 10^ ms~^ and 
9.1 X 10^ ms~^ calculated for the band of bismuth 
[25[. Our data are consistent with the extremely small 
effective mass of 0.002me (where rUe is the electron mass) 
observed in magneto-reflection measurements on samples 
with X = 11% The Dirac point in graphene, co- 

incidentally, has a band velocity {\vf\ ~ 10^ ms~^) [3^ 
comparable to what we observe for Bio.gSbo.i, but its 
spin-orbit coupling is several orders of magnitude weaker, 
and the only known method of inducing a gap in the 
Dirac spectrum of graphene is by coupling to an exter- 
nal chemical substrate [31]. The Bii-^^Sba^ series thus 
provides a rare opportunity to study relativistic Dirac 
Hamiltonian physics in a 3D condensed matter system 
where the intrinsic (rest) mass gap can be easily tuned. 

Studying the band dispersion perpendicular to the 
sample surface provides a way to differentiate bulk states 
from surface states in a 3D material. To visualize the 
neai-EF dispersion along the 3D L-X cut (X is a point 
that is displaced from L by a /c^ distance of 37r/c, where c 
is the lattice constant), in Fig. 2a we plot energy distribu- 
tion curves (EDCs), taken such that electrons at Ep have 
fixed in-plane momentum (kx^ky) — (L^^, Ly) = (0.8 A~^, 
0.0 A~^), as a function of photon energy {hu). There are 
three prominent features in the EDCs: a non-dispersing, 
kz independent, peak centered just below Ef dit about 
—0.02 eV; a broad non-dispersing hump centered near 
—0.3 eV; and a strongly dispersing hump that coincides 
with the latter near hu = 29 eV. To understand which 
bands these features originate from, we show ARPES in- 
tensity maps along an in-plane cut KMK (parallel to the 
ky direction) taken using hv values of 22 eV, 29 eV and 
35 eV, which correspond to approximate kz values of — 
0.3 A-\ Lz, and + 0.3 A"^ respectively (Fig.2b). At 
hv = 29 eV, the low energy ARPES spectral weight re- 
veals a clear A-shaped band close to Ep- As the photon 
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FIG. 2: The first 3D topological insulator (2007): Dirac-like dispersion signalling band inversion via spin- 
orbit interaction a, Schematic of an ARPES experimental geometry. The kinetic energy of photoelectrons and their angles 
of emission {0,(j)) determine its electronic structure, b, Energy and momentum analysis take place through a hemispherical 
analyzer and the spin analysis is performed using a Mott detector. Selected ARPES intensity maps of Bio.gSbo.i are shown 
along three /c-space cuts through the L point of the bulk 3D Brillouin zone (BZ). The presented data are taken in the third 
BZ with Lz = 2.9 with a photon energy of 29 eV. The cuts are along c, the ky direction, d, a direction rotated by 
approximately 10° from the ky direction, and e, the kx direction. Each cut shows a A-shaped bulk band whose tip lies below 
the Fermi level signalling a bulk gap. The (topological) surface states are denoted (T)SS and are all identified in Fig|3] [for 
further identification via theoretical calculations see Supplementary Materials (SM)]. f. Momentum distribution curves (MDCs) 
corresponding to the intensity map in c. h. Log scale plot of the MDCs corresponding to the intensity map in e. The red lines 
are guides to the eye for the bulk features in the MDCs. g. Schematic of the bulk 3D BZ of Bii-^Sb^ and the 2D BZ of the 
projected (111) surface. The high symmetry points F, M and K of the surface BZ are labeled. Schematic evolution of bulk 
band energies as a function of x is shown. The L band inversion transition occurs at x ~ 0.04, where a 3D gapless Dirac point 
is realized, and the composition we study here (for which x = 0.1) is indicated by the green arrow. A more detailed phase 
diagram based on our experiments is shown in FigU):. [Adapted from D. Hsieh et al., Nature 452, 970 (2008) (Completed and 
Submitted in 2007) ^14]]. 



energy is either increased or decreased from 29 eV, this 
intensity shifts to higher binding energies as the spectral 
weight evolves from the A-shaped into a U-shaped band. 
Therefore the dispersive peak in Fig. 2a comes from the 
bulk valence band, and for hu = 29 eV the high sym- 
metry point L = (0.8, 0, 2.9) appears in the third bulk 
BZ. In the maps of Fig. 2b with respective hu values of 
22 eV and 35 eV, overall weak features near Ep that 
vary in intensity remain even as the bulk valence band 
moves far below Ep. The survival of these weak features 
over a large photon energy range (17 to 55 eV) supports 
their surface origin. The non-dispersing feature centered 
near —0.3 eV in Fig. 2a comes from the higher binding 



energy (valence band) part of the full spectrum of sur- 
face states, and the weak non-dispersing peak at —0.02 
eV reflects the low energy part of the surface states that 
cross Ef away from the M point and forms the surface 
Fermi surface (Fig. 2c). 
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FIG. 3: The first 3D topological insulator (2007): Topological Surface States and electronic band dispersion 
along the kz-direction in momentum space. Surface states are experimentally identified by studying their out-of-plane 
momentum dispersion through the systematic variation of incident photon energy, a, Energy distribution curves (EDCs) of 
Bio.gSbo.i with electrons at the Fermi level (Ep) maintained at a fixed in-plane momentum of {kx=0.8 A~^, ky^O.O A~^) are 
obtained as a function of incident photon energy to identify states that exhibit no dispersion perpendicular to the (lll)-plane 
along the direction shown by the double-headed arrow labeled "3" in the inset (see SM). Selected EDO data sets with photon 
energies of 28 eV to 32 eV in steps of 0.5 eV are shown for clarity. The non-energy dispersive (kz independent) peaks near Ef 
are the topological surface states (TSS). b, ARPES intensity maps along cuts parallel to ky taken with electrons at Ef fixed 
at kx — 0.8 with respective photon energies of hv — 22 eV, 29 eV and 35 eV (for a conversion map from photon energy 
to kz see SM). The faint A-shaped band at hv = 22 eV and hv = 35 eV shows some overlap with the bulk valence band at L 
{hv — 29 eV), suggesting that it is a resonant surface state degenerate with the bulk state in some limited k-range near Ef- 
The flat band of intensity centered about —2 eV in the hf — 22 eV scan originates from Bi 5d core level emission from second 
order light, c, Projection of the bulk BZ (black lines) onto the (111) surface BZ (green lines). Overlay (enlarged in inset) 
shows the high resolution Fermi surface (FS) of the metallic SS mode, which was obtained by integrating the ARPES intensity 
(taken with hv — 20 eV) from —15 meV to 10 meV relative to Ef- The six tear-drop shaped lobes of the surface FS close to F 
(center of BZ) show some intensity variation between them that is due to the relative orientation between the axes of the lobes 
and the axis of the detector slit. The six- fold symmetry was however confirmed by rotating the sample in the kx — ky plane. 
EDCs corresponding to the cuts A, B and C are also shown; these confirm the gap less character of the surface states in bulk 
insulating Bio.gSbo.i. [Adapted from D. Hsieh et al., Nature 452, 970 (2008) (Completed and Submitted in 2007) [14]]. 



WINDING NUMBER COUNT: COUNTING OF 
SURFACE FERMI SURFACES ENCLOSING 
KRAMERS POINTS TO IDENTIFY 
TOPOLOGICALLY NON-TRIVIAL SURFACE 
SPIN-TEXTURED STATES 



Having established the existence of an energy gap in 
the bulk state of Bio.gSbo.i (Figs 1 and 2) and observed 



linearly dispersive bulk bands uniquely consistent with 
strong spin-orbit coupling model calculations [25*428| (see 
SM for full comparison with theoretical calculation), we 
now discuss the topological character of its surface states, 
which are found to be gapless (Fig. 2c). In general, the 
states at the surface of spin-orbit coupled compounds 
are allowed to be spin split owing to the loss of space 
inversion symmetry [£^(/c,t) = ^^t)]- However, as 
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FIG. 4: The first 3D topological insulator (2007): The topological gapless surface states in bulk insulating 
(observed via the bulk band dispersion gap) Bio.gSbo.i. a, The surface band dispersion second derivative image (SDI) 
of Bio.gSbo.i along T — M. The shaded white area shows the projection of the bulk bands based on ARPES data, as well as 
a rigid shift of the tight binding bands to sketch the unoccupied bands above the Fermi level. To maintain high momentum 
resolution, data were collected in two segments of momentum space, then the intensities were normalized using background 
level above the Fermi level. A non-intrinsic flat band of intensity near Ef generated by the SDI analysis was rejected to isolate 
the intrinsic dispersion. The Fermi crossings of the surface state are denoted by yellow circles, with the band near —kx ~ 0.5 
counted twice owing to double degeneracy. The red lines are guides to the eye. An in-plane rotation of the sample by 60° 
produced the same surface state dispersion. The EDCs along F — M are shown to the right. There are a total of five crossings 
from F — M which indicates that these surface states are topologically non-trivial. The number of surface state crossings in a 
material (with an odd number of Dirac points) is related to the topological Z2 invariant (see text), b. The resistivity curves 
of Bi and Bio.gSbo.i reflect the contrasting transport behaviours. The presented resistivity curve for pure bismuth has been 
multiplied by a factor of 80 for clarity, c. Schematic variation of bulk band energies of Bii-^cSbcc as a function of x (based on 
band calculations and on 10, 24]). Bio.gSbo.i is a direct gap bulk Dirac point insulator well inside the inverted-band regime, 
and its surface forms a "topological metal" - the 2D analogue of the ID edge states in quantum spin Hall systems, d, ARPES 
intensity integrated within ±10 meV of originating solely from the surface state crossings. The image was plotted by 
stacking along the negative kx direction a series of scans taken parallel to the ky direction, e. Outline of Bio.gSbo.i surface 
state ARPES intensity near Ep measured in d. White lines show scan directions "1" and "2". f. Surface band dispersion 
along direction "1" taken with hv = 28 eV and the corresponding EDCs (g). The surface Kramers degenerate point, critical 
in determining the topological Z2 class of a band insulator, is clearly seen at M, approximately 15 ± 5 meV below Ep. (We 
note that the scans are taken along the negative kx direction, away from the bulk L point.) h. Surface band dispersion along 
direction "2" taken with hu = 28 eV and the corresponding EDCs (i). This scan no longer passes through the M-point, and 
the observation of two well separated bands indicates the absence of Kramers degeneracy as expected, which cross-checks the 
resuh in (a). [Adapted from D. Hsieh et a/.. Nature 452, 970 (2008) (Completed and Submitted in 2007)[iJ 



required by Kramers' theorem, this splitting must go to 
zero at the four time reversal invariant momenta (TRIM) 
in the 2D surface BZ. As discussed in [10, 12], along a 
path connecting two TRIM in the same BZ, the Fermi 
energy inside the bulk gap will intersect these singly de- 
generate surface states either an even or odd number of 



times. When there are an even number of surface state 
crossings, the surface states are topologically trivial be- 
cause weak disorder (as may arise through alloying) or 
correlations can remove pairs of such crossings by push- 
ing the surface bands entirely above or below Ef. When 
there are an odd number of crossings, however, at least 
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FIG. 5: Spin texture of a topological insulator encodes Z2 topological order of the bulk (A) Schematic sketches of 
the bulk Brillouin zone (BZ) and (111) surface BZ of the Bii-^Sbx crystal series. The high symmetry points (L,H,T,r,r,M,K) 
are identified. (B) Schematic of Fermi surface pockets formed by the surface states (SS) of a topological insulator that carries 
a Berry's phase. (C) Partner switching band structure topology: Schematic of spin-polarized SS dispersion and connectivity 
between F and M required to realize the FS pockets shown in panel- (B). La and Ls label bulk states at L that are antisymmetric 
and symmetric respectively under a parity transformation (see text). (D) Spin-integrated ARPES intensity map of the SS of 
Bio.9iSbo.o9 at Ef- Arrows point in the measured direction of the spin. (E) High resolution ARPES intensity map of the SS 
at that enclose the Mi and M2 points. Corresponding band dispersion (second derivative images) are shown below. The 
left right asymmetry of the band dispersions are due to the slight offset of the alignment from the f-Mi(M2) direction. (F) 
Surface band dispersion image along the F-M direction showing five Fermi level crossings. The intensity of bands 4,5 is scaled 
up for clarity (the dashed white lines are guides to the eye). The schematic projection of the bulk valence and conduction 
bands are shown in shaded blue and purple areas. (G) Spin-resolved momentum distribution curves presented at Eb = —25 
meV showing single spin degeneracy of bands at 1, 2 and 3. Spin up and down correspond to spin pointing along the +y and 
-y direction respectively. (H) Schematic of the spin-polarized surface FS observed in our experiments. It is consistent with a 
— 1 topology (compare (B) and (H)). [Adapted from D. Hsieh et al., Science 323, 919 (2009) (Completed and Submitted 
in 2008) [ll]. 



one surface state must remain gapless, which makes it 
non-trivial [lo|, [HI, Hsj . The existence of such topologi- 
cally non-trivial surface states can be theoretically pre- 
dicted on the basis of the bulk band structure only, us- 
ing the Z2 invariant that is related to the quantum Hall 
Chern number [8[. Materials with band structures with 
Z2 = +1 are ordinary Bloch band insulators that are 
topologically equivalent to the filled shell atomic insula- 
tor, and are predicted to exhibit an even number (includ- 
ing zero) of surface state crossings. Materials with bulk 
band structures with Z2 = — 1 on the other hand, which 
are expected to exist in rare systems with strong spin- 
orbit coupling acting as an internal quantizing magnetic 
field on the electron system [32|, and inverted bands at 
an odd number of high symmetry points in their bulk 3D 
BZs, are predicted to exhibit an odd number of surface 
state crossings, precluding their adiabatic continuation to 
the atomic insulator [9- 13, "23]. Such "topological quan- 
tum Hall metals" JJ] cannot be realized in a purely 



2D electron gas system such as the one realized at the 
interface of GaAs/GaAlAs systems. 

In our experimental case, namely the (111) surface of 
Bio.gSbo.i, the four TRIM are located at f and three M 
points that are rotated by 60° relative to one another. 
Owing to the three-fold crystal symmetry (A7 bulk struc- 
ture) and the observed mirror symmetry of the surface 
Fermi surface across kx = (Fig. 2), these three M points 
are equivalent (and we henceforth refer to them as a sin- 
gle point, M). The mirror symmetry [E{ky) = E{—ky)] 
is also expected from time reversal invariance exhibited 
by the system. The complete details of the surface state 
dispersion observed in our experiments along a path con- 
necting f and M are shown in Fig. 3a; finding this infor- 
mation is made possible by our experimental separation 
of surface states from bulk states. As for bismuth (Bi), 
two surface bands emerge from the bulk band continuum 
near F to form a central electron pocket and an adjacent 
hole lobe (33l-[35|. It has been established that these two 
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FIG. 6: Spin-momentum locked helical fermions reveal 
topological-order of the bulk: Spin-momentum locked he- 
lical surface Dirac fermions are hallmark signatures of topo- 
logical insulators, (a) ARPES data for Bi2Se3 reveals surface 
electronic states with a single spin-polarized Dirac cone. The 
surface Fermi surface (b) exhibits a chiral left-handed spin 
texture. Data is from a sample doped to a Fermi level [thin 
line in (a)] that is in the bulk gap. (c) Surface electronic struc- 
ture of Bi2Se3 computed in the local density approximation. 
The shaded regions describe bulk states, and the red lines 
are surface states, (d) Schematic of the spin polarized surface 
state dispersion in the 612X3 topological insulators. [Adapted 
from Y. Xia et al., Nature Phys. 5, 398 (2009). (Completed 
and Submitted in 2008) 15]]. 



bands result from the spin-splitting of a surface state 
and are thus singly degenerate [35|, [36|. On the other 
hand, the surface band that crosses Ep at —kx ~ 0.5 
A~^, and forms the narrow electron pocket around M, 
is clearly doubly degenerate, as far as we can determine 
within our experimental resolution. This is indicated by 
its splitting below Ep between —kx ~ 0.55 and M, 
as well as the fact that this splitting goes to zero at M 
in accordance with Kramers theorem. In semimetallic 
single crystal bismuth, only a single surface band is ob- 
served to form the electron pocket around M [s^, [ssf . 
Moreover, this surface state overlaps, hence becomes de- 
generate with, the bulk conduction band at L (L projects 
to the surface M point) owing to the semimetallic char- 
acter of Bi (Fig. 3b). In Bio.gSbo.i on the other hand, the 
states near M fall completely inside the bulk energy gap 
preserving their purely surface character at M (Fig. 3a). 
The surface Kramers doublet point can thus be defined 
in the bulk insulator (unlike in Bi [33l438j ) and is ex- 
perimentally located in Bio.gSbo.i samples to lie approx- 
imately 15 ± 5 meV below Ep at k = M (Fig. 3a). For 
the precise location of this Kramers point, it is impor- 
tant to demonstrate that our alignment is strictly along 



the f — M line. To do so, we contrast high resolution 
ARPES measurements taken along the T — M line with 
those that are slightly offset from it (Fig.3e). Figs 3f-i 
show that with ky offset from the Kramers point at M by 
less than 0.02 A~^, the degeneracy is lifted and only one 
band crosses Ej? to form part of the bow-shaped elec- 
tron distribution (Fig. 3d). Our finding of five surface 
state crossings (an odd rather than an even number) be- 
tween f and M (Fig. 3a), confirmed by our observation 
of the Kramers degenerate point at the TRIM, indicates 
that these gapless surface states are topologically non- 
trivial. This corroborates our bulk electronic structure 
result that Bio.gSbo.i is in the insulating band- inverted 
{Z2 = — 1) regime (Fig. 3c), which contains an odd num- 
ber of bulk (gapped) Dirac points, and is topologically 
analogous to an integer quantum spin Hall insulator. 

Our experimental results taken collectively strongly 
suggest that Bio.gSbo.i is quite distinct from graphene 
[30, ^9] and represents a novel state of quantum mat- 
ter: a strongly spin-orbit coupled insulator with an odd 
number of Dirac points with a negative Z2 topological 
Hall phase, which realizes the "parity anomaly without 
Fermion doubling". Our work further demonstrates a 
general methodology for possible future investigations of 
novel topological orders in exotic quantum matter. 



SPIN-RESOLVING THE SURFACE STATES TO 
IDENTIFY THE NON-TRIVIAL TOPOLOGICAL 
PHASE AND ESTABLISH A 2D HELICAL 
METAL PROTECTED FROM 
BACKSCATTERING 

Strong topological materials are distinguished from or- 
dinary materials such as gold by a topological quantum 
number, i^o = 1 or respectively @, [n|. For Bii-a^Sb^^, 
theory has shown that uq is determined solely by the char- 
acter of the bulk electronic wave functions at the L point 
in the three-dimensional (3D) Brillouin zone (BZ). When 
the lowest energy conduction band state is composed of 
an antisymmetric combination of atomic p-type orbitals 
{La) and the highest energy valence band state is com- 
posed of a symmetric combination (Ls), then i^o = 1, and 
vice versa for i^o = [10]. Although the bonding nature 
(parity) of the states at L is not revealed in a measure- 
ment of the bulk band structure, the value of can be 
determined from the spin-textures of the surface bands 
that form when the bulk is terminated. In particular, a 
i^o = 1 topology requires the terminated surface to have 
a Fermi surface (FS) that supports a non-zero Berry's 
phase (odd as opposed to even multiple of tt), which is 
not realizable in an ordinary spin-orbit material. 

In a general inversion symmetric spin-orbit insulator, 
the bulk states are spin degenerate because of a combi- 
nation of space inversion symmetry [E{k,^) = E{—k,^)] 
and time reversal symmetry [E{k^'\') = E{—k^l.)]. Be- 
cause space inversion symmetry is broken at the termi- 
nated surface, the spin degeneracy of surface bands can 
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FIG. 7: Topological character of parent compound revealed on the (111) surface states. Schematic of the bulk band 
structure (shaded areas) and surface band structure (red and blue lines) of Sb near Ep for a (A) topologically non-trivial and 
(B) topological trivial (gold-like) case, together with their corresponding surface Fermi surfaces are shown. (C) Spin-integrated 
ARPES spectrum of Sb(lll) along the F-M direction. The surface states are denoted by SS, bulk states by BS, and the hole-like 
resonance states and electron- like resonance states by h RS and e~ RS respectively. (D) Calculated surface state band structure 
of Sb(lll) based on the methods in [li,!!^. The continuum bulk energy bands are represented with pink shaded regions, and 
the lines show the discrete bands of a 100 layer^slab. The red and blue single bands, denoted Si and 112, are the surface states 
bands with spin polarization (P) oc -\-y and (P) oc —y respectively. (E) ARPES intensity map of Sb(lll) at Ef in the kx-ky 
plane. The only one FS encircling F seen in the data is formed by the inner V-shaped SS band seen in panel- (C) and (F). 
The outer V-shaped band bends back towards the bulk band best seen in data in panel- (F). (F) ARPES spectrum of Sb(lll) 
along the F-K direction shows that the outer V-shaped SS band merges with the bulk band. (G) Schematic of the surface FS 
of Sb(lll) showing the pockets formed by the surface states (unfilled) and the resonant states (blue and purple). The purely 
surface state Fermi pocket encloses only one Kramers degenerate point (^t), namely, f (=^t), therefore consistent with the 
— 1 topological classification of Sb which is different from Au (compare (B) and (G)). As discussed in the text, the hRS and 
e~RS count trivially. [Adapted from D. Hsieh et a/.. Science 323, 919 (2009) (Completed and Submitted in 2008) O]. 



be lifted by the spin-orbit interaction [19-21]. However, 
according to Kramers theorem [16], they must remain 
spin degenerate at four special time reversal invariant 
momenta (kr = T, M) in the surface BZ [11], which for 
the (111) surface of Bii_a^Sba^ are located at f and three 
equivalent M points [see FigEfA)]. 

Depending on whether equals or 1, the Fermi sur- 
face pockets formed by the surface bands will enclose the 
four kT an even or odd number of times respectively. If 
a Fermi surface pocket does not enclose kr (= f , M), it 
is irrelevant for the topology [lo|, I1q| . Because the wave 
function of a single electron spin acquires a geometric 
phase factor of tt (isl ? ] as it evolves by 360° in mo- 
mentum space along a Fermi contour enclosing a A:t, an 
odd number of Fermi pockets enclosing kx in total im- 
plies a TT geometrical (Berry's) phase [10]. In order to 
realize a tt Berry's phase the surface bands must be spin- 
polarized and exhibit a partner switching [10] dispersion 
behavior between a pair of kr- This means that any pair 
of spin-polarized surface bands that are degenerate at F 
must not re-connect at M, or must separately connect 
to the bulk valence and conduction band in between f 
and M. The partner switching behavior is realized in Fig. 



[5][C) because the spin down band connects to and is de- 
generate with different spin up bands at F and M. The 
partner switching behavior is realized in Fig. [Tl^A) be- 
cause the spin up and spin down bands emerging from 
f separately merge into the bulk valence and conduction 
bands respectively between f and M. 

We first investigate the spin properties of the topo- 
logical insulator phase [l6|. Spin-integrated ARPES 5l| 
intensity maps of the (111) surface states of insulating 
Bii_^Sb^ taken at the Fermi level {Ep) [Figs[5];D)&(E)] 
show that a hexagonal FS encloses F, while dumbbell 
shaped FS pockets that are much weaker in intensity en- 
close M. By examining the surface band dispersion be- 
low the Fermi level [FiglSfF)] it is clear that the central 
hexagonal FS is formed by a single band (Fermi crossing 
1) whereas the dumbbell shaped FSs are formed by the 
merger of two bands (Fermi crossings 4 and 5) p^ . 

This band dispersion resembles the partner switching 
dispersion behavior characteristic of topological insula- 
tors. To check this scenario and determine the topo- 
logical index z^o, we have carried out spin-resolved pho- 
toemission spectroscopy. FigEfG) shows a spin-resolved 
momentum distribution curve taken along the f -M direc- 
tion at a binding energy Eb = —25 meV [FigEfG)]. The 
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FIG. 8: Observation of bulk band inversion across a 
topological quantum phase transition. Energy levels of 
r4 . (blue circles) and (green squares) bands are obtained 
from ARPES measurements as a function of composition S. 
CB: conduction band; VB: valence band. Parity eigenvalues 
(+ or ) of Bloch states are shown. The topological invariants, 
Vi, obtained from the parity eigenvalues are presented as [O/tt 
= 1^0 ;iy 1,1^2, 1^3] where = tvuo is the axion angle [44] and is 
the strong invariant. [Adapted from S.-Y. Xu et al., Science 
332 560 (2011). [50(] 



data reveal a clear difference between the spin-up and 
spin-down intensities of bands 1, 2 and 3, and show that 
bands 1 and 2 have opposite spin whereas bands 2 and 
3 have the same spin (detailed analysis discussed later 
in text). The former observation confirms that bands 1 
and 2 form a spin-orbit split pair, and the latter obser- 
vation suggests that bands 2 and 3 (as opposed to bands 
1 and 3) are connected above the Fermi level and form 
one band. This is further confirmed by directly imag- 
ing the bands through raising the chemical potential via 
doping (SM). Irrelevance of bands 2 and 3 to the topol- 
ogy is consistent with the fact that the Fermi surface 
pocket they form does not enclose any kr- Because of 
a dramatic intrinsic weakening of signal intensity near 
crossings 4 and 5, and the small energy and momentum 
splitting of bands 4 and 5 lying at the resolution limit 
of modern spin-resolved ARPES spectrometers, no con- 
clusive spin information about these two bands can be 
drawn from the methods employed in obtaining the data 
sets in Figs [5jG)&(H). However, whether bands 4 and 
5 are both singly or doubly degenerate does not change 
the fact that an odd number of spin-polarized FSs en- 
close the /ct, which provides evidence that Bii_a^Sba^ has 
1^0 = 1 and that its surface supports a non-trivial Berry's 
phase. This directly implies an absence of backscattering 
in electronic transport along the surface (FigITT]). 

Shortly after the discovery of Bii-a^Sb^^, physicists 



sought to find a simpler version of a 3D topological insula- 
tor consisting of a single surface state instead of five. This 
is because the surface structure of Bii_a;Sba; was rather 
complicated and the band gap was rather small. This mo- 
tivated a search for topological insulators with a larger 
band gap and simpler surface spectrum. A second gen- 
eration of 3D topological insulator materials, especially 
Bi2Se3, offers the potential for topologically protected be- 
havior in ordinary crystals at room temperature and zero 
magnetic field. Starting in 2008, work by the Princeton 
group used spin- ARPES and first-principles calculations 
to study the surface band structure of Bi2Se3 and observe 
the characteristic signature of a topological insulator in 
the form of a single Dirac cone that is spin-polarized (Figs 
[6] and ITQI) such that it also carries a non-trivial Berry's 
phase jlBI. fl7|. Concurrent theoretical work by [19] used 
electronic structure methods to show that Bi2Se3 is just 
one of several new large band-gap topological insulators. 
These other materials were soon after also identified us- 
ing this ARPES technique we describe 0,0. 

The Bi2Se3 surface state is found from spin- ARPES 
and theory to be a nearly idealized single Dirac cone as 
seen from the experimental data in Figsl6] and [lH An 
added advantage is that Bi2Se3 is stoichiometric (i.e., a 
pure compound rather than an alloy such as Bii-a^Sb^;) 
and hence can be prepared, in principle, at higher purity. 
While the topological insulator phase is predicted to be 
quite robust to disorder, many experimental probes of the 
phase, including ARPES of the surface band structure, 
are clearer in high-purity samples. Finally and perhaps 
most important for applications, Bi2Se3 has a large band 
gap of around 0.3 eV (3600 K). This indicates that in its 
high-purity form Bi2Se3 can exhibit topological insulator 
behavior at room temperature and greatly increases the 
potential for applications, which we discuss in greater 
depth later in the review. 



IDENTIFYING THE ORIGIN OF 3D 
TOPOLOGICAL ORDER VIA A BULK BAND 
GAP INVERSION TRANSITION 

We investigated the quantum origin of topological or- 
der in this class of materials. It has been theoretically 
speculated that the novel topological order originates 
from the parities of the electrons in pure Sb and not 
Bi [io|, 111. It was also noted that the origin of 
the topological effects can only be tested by measuring 
the spin-texture of the Sb surface, which has not been 
measured. Based on quantum oscillation and magneto- 
optical studies, the bulk band structure of Sb is known 
to evolve from that of insulating Bii-a^Sb^^ through the 
hole-like band at H rising above^i? and the electron-like 
band at L sinking below Ef j2^. The relative energy 
ordering of the La and Ls states in Sb again determines 
whether the surface state pair emerging from f switches 
partners [Fig|7l^A)] or not [Fig|7l^B)] between f and M, 
and in turn determines whether they support a non-zero 
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FIG. 9: Spin-texture of topological surface states and topological chirality. (A) Experimental geometry of the spin- 
resolved ARPES study. At normal emission {0 = 0°), the sensitive y'-axis of the Mott detector is rotated by 45° from the 
sample f to — M (|| —x) direction, and the sensitive z^-axis of the Mott detector is parallel to the sample normal (|| z). (B) 
Spin-integrated ARPES spectrum of Sb(lll) along the — M-F-M direction. The momentum splitting between the band minima 
is indicated by the black bar and is approximately 0.03 A~^. A schematic of the spin chirality of the central FS based on 
the spin- resolved ARPES results is shown on the right. (C) Momentum distribution curve of the spin averaged spectrum at 
Eb = —30 meV (shown in (B) by white line), together with the Lorentzian peaks of the fit. (D) Measured spin polarization 
curves (symbols) for the detector y' and z' components together with the fitted lines using the two-step fitting routine [4^. 
(E) Spin-resolved spectra for the sample y component based on the fitted spin polarization curves shown in (D). Up (down) 
triangles represent a spin direction along the -\-{-)y direction. (F) The in-plane and out-of-plane spin polarization components 
in the sample coordinate frame obtained from the spin polarization fit. Overall spin-resolved data and the fact that the surface 
band that forms the central electron pocket has (P) oc —y along the -\-kx direction, as in (E), suggest a left-handed chirality 
(schematic in (B) and see text for details). [Adapted from D. Hsieh et al., Science 323, 919 (2009) ^16.]]. 



Berry's phase. 

In a conventional spin-orbit metal such as gold, a free- 
electron like surface state is split into two parabolic spin- 
polarized sub-bands that are shifted in ^-space relative 
to each other [46j. Two concentric spin-polarized Fermi 
surfaces are created, one having an opposite sense of in- 
plane spin rotation from the other, that enclose f . Such 
a Fermi surface arrangement, like the schematic shown 
in figure [7JB), does not support a non-zero Berry's phase 
because the kr are enclosed an even number of times (2 
for most known materials). 

However, for Sb, this is not the case. Figure [H^C) 
shows a spin-integrated ARPES intensity map of Sb(lll) 
from F to M. By performing a systematic incident pho- 
ton energy dependence study of such spectra, previously 
unavailable with He lamp sources [47], it is possible to 
identify two V-shaped surface states (SS) centered at 
F, a bulk state located near ka: = -0.25 A-1 and res- 
onance states centered about kx = 0.25 A~^ and M that 
are hybrid states formed by surface and bulk states 



(SM). An examination of the ARPES intensity map of 
the Sb(lll) surface and resonance states at Ef [FiglTjE)] 
reveals that the central surface FS enclosing f is formed 
by the inner V-shaped SS only. The outer V-shaped SS 
on the other hand forms part of a tear-drop shaped FS 
that does not enclose f , unlike the case in gold. This 
tear-drop shaped FS is formed partly by the outer V- 
shaped SS and partly by the hole-like resonance state. 
The electron-like resonance state FS enclosing M does 
not affect the determination of uq because it must be dou- 
bly spin degenerate (SM). Such a FS geometry [FigEjG)] 
suggests that the V-shaped SS pair may undergo a part- 
ner switching behavior expected in FigCjA). This behav- 
ior is most clearly seen in a cut taken along the f -K direc- 
tion since the top of the bulk valence band is well below 
Ef [Figlll^F)] showing only the inner V-shaped SS cross- 
ing Ef while the outer V-shaped SS bends back towards 
the bulk valence band near kx = 0.1 A~^ before reach- 
ing Ef' The additional support for this band dispersion 
behavior comes from tight binding surface calculations 
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FIG. 10: First detection of Topological- Order: spin-momentum locking of spin-helical Dirac electrons in Bi2Se3 
and Bi2Te3 using spin-resolved ARPES. (a) ARPES intensity map at Ep of the (111) surface of tuned Bi2-(5Ca5Se3 
(see text) and (b) the (111) surface of Bi2Te3. Red arrows denote the direction of spin around the Fermi surface, (c) ARPES 
dispersion of tuned Bi2-5Ca5Se3 and (d) Bi2Te3 along the kx cut. The dotted red hues are guides to the eye. The shaded regions 
in (c) and (d) are our calculated projections of the bulk bands of pure Bi2Se3 and Bi2Te3 onto the (111) surface respectively, 
(e) Measured y component of spin-polarization along the F-M direction at Eb = -20 meV, which only cuts through the surface 
states. Inset shows a schematic of the cut direction, (f) Measured x (red triangles) and z (black circles) components of 
spin-polarization along the F-M direction at Eb — -20 meV. Error bars in (e) and (f) denote the standard deviation of Px,y,z, 
where typical detector counts reach 5 x 10^; Solid lines are numerical fits [43]. (g) Spin-resolved spectra obtained from the 
y component spin polarization data. The non-Lorentzian lineshape of the ly and ly curves and their non-exact merger at 
large \kx\ is due to the time evolution of the surface band dispersion, which is the dominant source of statistical uncertainty, 
a.u., arbitrary units, (h) Fitted values of the spin polarization vector P = (Sx^Sy^Sz) are (sin(90°)cos(-95°), sin(90°)sin(-95°), 
cos(90°)) for electrons with -\-kx and (sin(86°)cos(85°), sin(86°)sin(85°), cos(86°)) for electrons with -kx, which demonstrates 
the topological helicity of the spin-Dirac cone. The angular uncertainties are of order ±10° and the magnitude uncertainty is 
of order ±0.15. [Adapted from D. Hsieh et al., Nature 460, 1101 (2009). [17]]. 



on Sb [FiglTl^D)], which closely match with experimental 
data below Ep. Our observation of a single surface band 
forming a FS enclosing F suggests that pure Sb is likely 
described by i^o = 1, and that its surface may support a 
Berry's phase. 

Confirmation of a surface tt Berry's phase rests criti- 
cally on a measurement of the relative spin orientations 
(up or down) of the SS bands near f so that the partner 
switching is indeed realized, which cannot be done with- 
out spin resolution. Spin resolution was achieved using 
a Mott polarimeter that measures two orthogonal spin 
components of a photoemitted electron [48, 49]. These 
two components are along the and directions of 
the Mott coordinate frame, which lie predominantly in 
and out of the sample (111) plane respectively. Each of 
these two directions represents a normal to a scattering 
plane defined by the photoelectron incidence direction 



on a gold foil and two electron detectors mounted on 
either side (left and right) [FiglQl^A)]. Strong spin-orbit 
coupling of atomic gold is known to create an asymme- 
try in the scattering of a photoelectron off the gold foil 
that depends on its spin component normal to the scat- 
tering plane [49|. This leads to an asymmetry between 
the left intensity {ly, and right intensity (/^ given 
by A^.,,. = (7^^,^,, - + ly'^z')^ which is re- 

lated to the spin polarization Py' ^z' = i^/Seff) x ^y'.z' 
through the Sherman function Seff = 0.085 [48|, |49|. 
Spin-resolved momentum distribution curve data sets of 
the SS bands along the -M-f-M cut at Eb = -30 meV 
[FiglD^B)] are shown for maximal intensity. Figure [9l^D) 
displays both y' and z' polarization components along 
this cut, showing clear evidence that the bands are spin 
polarized, with spins pointing largely in the (111) plane. 
In order to estimate the full 3D spin polarization vec- 
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FIG. 11: Helical spin texture naturally leads to absence of elastic backscattering for surface transport: No "U" 
turn on a 3D topological insulator surface, (a) Our measurement of a helical spin texture in both Bii-^Sbx and in 
Bi2Se3 directly shows that there is (b) an absence of backscattering. (c) ARPES measured FSs are shown with spin directions 
based on polarization measurements. L(R)HC stands for left(right)-handed chirality. (d) Spin independent and spin dependent 
scattering profiles on FSs in (c) relevant for surface quasi-particle transport are shown which is sampled by the quasi-particle 
interference (QPI) modes. [Adapted from S.-Y. Xu et a/., Science 332 560 (2011). [SC 



tors from a two component measurement (which is not 
required to prove the partner switching or the Berry's 
phase), we fit a model polarization curve to our data fol- 
lowing the recent demonstration in Ref-[45], which takes 
the polarization directions associated with each momen- 
tum distribution curve peak [FiglHl^C)] as input param- 
eters, with the constraint that each polarization vector 
has length one (in angular momentum units of fi/2). Our 
fitted polarization vectors are displayed in the sample 
(x^y^z) coordinate frame [Figl9fF)], from which we de- 
rive the spin-resolved momentum distribution curves for 
the spin components parallel (/^) and anti-parallel (/^) 
to the y direction (SM) as shown in figure [91(E). There is 
a clear difference in and at each of the four momen- 
tum distribution curve peaks indicating that the surface 
state bands are spin polarized [FiglQl^E)], which is possi- 
ble to conclude even without a full 3D fitting. Each of the 
pairs /2//1 and rl/r2 have opposite spin, consistent with 
the behavior of a spin split pair, and the spin polariza- 
tion of these bands are reversed on either side of f in ac- 



cordance with the system being time reversal symmetric 
[E{k,^) = E{-k,i)] [FigijF)]. The measured spin tex- 
ture of the Sb(lll) surface states (Fig|9]), together with 
the connectivity of the surface bands (Fig|7|), uniquely 
determines its belonging to the uq = 1 class. Therefore 
the surface of Sb carries a non-zero (tt) Berry's phase via 
the inner V-shaped band and pure Sb can be regarded 
as the parent metal of the Bii-^^Sb^; topological insu- 
lator class, in other words, the topological order origi- 
nates from the Sb wave functions. A recent work has 
demonstrated a topological quantum phase transition as 
a function of chemical composition from a non-inverted 
to an inverted semiconductor as a clear example of the 
origin of topological order (FigHV|). 

Our spin polarized measurement methods (Fig|5] and 
[9|) uncover a new type of topological quantum number 
riM which provides information about the chirality prop- 
erties. Topological band theory suggests that the bulk 
electronic states in the mirror {ky = 0) plane can be 
classified in terms of a number um (=il) that describes 
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3D Vector spin texture 




Momentum 



FIG. 12: Spin texture evolution of topological surface 
bands as a function of energy away from the Dirac 
node (left axis) and geometrical warping factor uj (right axis). 
The warping factor is defined as cj = ^^(F 



where uj 

2 — v3 



kF{T-M)-kF{T-K) 
-M) + kF{T-K) 



0, cj = 1, and cj > 1 implies circular, hexag- 
onal and snowflake-shaped FSs respectively. The sign of uj in- 
dicates texture chirality for LHC (+) or RHC (-). Insets: out- 
of-plane 3D spin-polarization measurements at corresponding 
FSs. [Adapted from S.-Y. Xu et a/., Science 332 560 (2011). 
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FIG. 13: Robustness against disorder Spin split surface 
states survive alloying disorder in Sbo.gBio.i. (A) ARPES 
intensity map at Ef of single crystal Sbo.gBio.i (111) in the 
kx-ky plane taken using 20 eV photons. (B) ARPES intensity 
map of Sbo.gBio.i (111) along the F-M direction taken with hu 
— 22 eV photons. The band dispersion is not symmetric about 
F because of the three-fold rotational symmetry of the bulk 
states about the (111) axis. [Adapted from D. Hsieh et al., 
Science 323, 919 (2009) 0]. 



the handedness (either left or right handed) or chirality 
of the surface spins which can be directly measured or 
seen in spin-resolved experiments [40|. We now deter- 
mine the value of um from our data. From figure [5l it 
is seen that a single (one) surface band, which switches 
partners at M, connects the bulk valence and conduction 
bands, so \nM\ = 1 (SM). The sign of um is related to 
the direction of the spin polarization (P) of this band 
(ioj . which is constrained by mirror symmetry to point 
along dz^. Since the central electron-like FS enclosing f 
intersects six mirror invariant points [see FiglHl^B)], the 
sign of Um distinguishes two distinct types of handed- 
ness for this spin polarized FS. Figures [5fF) and [9] show 
that for both Bii-^^Sba^ and Sb, the surface band that 
forms this electron pocket has (P) oc —y along the kx 
direction, suggesting a left-handed rotation sense for the 
spins around this central FS thus um = — 1- Therefore, 
both insulating Bii-^^Sba^ and pure Sb possess equivalent 
chirality properties — a definite spin rotation sense (left- 
handed chirality, see Figl9l^B)) and a topological Berry's 
phase. Recently a chirality transition across the surface 
Dirac point of a 3D topological insulator has also been 
observed using spin- ARPES [50|] {Fig^. 

These spin-resolved experimental measurements reveal 
an intimate and straightforward connection between the 
topological numbers {uq^ um) and the physical observ- 
ables. The determines whether the surface electrons 
support a non-trivial Berry's phase, and if they do, the 
Um determines the spin handedness of the Fermi surface 
that manifests this Berry's phase. The 2D Berry's phase 
is a critical signature of topological order and is not real- 
izable in isolated 2D electron systems, nor on the surfaces 
of conventional spin-orbit or exchange coupled magnetic 
materials. A non-zero Berry's phase is known, theoret- 
ically, to protect an electron system against the almost 
universal weak-localization behavior in their low temper- 
ature transport 0, ^| and is expected to form the key 
element for fault-tolerant computation schemes 0, [52|, 
because the Berry's phase is a geometrical agent or mech- 
anism for protection against quantum decoherence [53^ . 
Its remarkable realization on the Bii-^^Sba^ surface repre- 
sents an unprecedented example of a 2D tt Berry's phase, 
and opens the possibility for building realistic prototype 
systems to test quantum computing modules. In gen- 
eral, our results demonstrate that spin- ARPES is a pow- 
erful probe of 3D topological order, which opens up a 
new search front for topological materials for novel spin- 
devices and fault-tolerant quantum computing. 



TOPOLOGICAL PROTECTION AND 
TUNABILITY OF THE SURFACE STATES OF A 
3D TOPOLOGICAL INSULATOR 

The predicted topological protection of the surface 
states of Sb implies that their metallicity cannot be de- 
stroyed by weak time reversal symmetric perturbations. 
In order to test the robustness of the measured gapless 
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FIG. 14: Observation of room temperature (300K) topological order without applied magnetic field in Bi2Se3: 

(a) Crystal momentum integrated ARPES data near Fermi level exhibit linear fall-off of density of states, which, combined with 
the spin-resolved nature of the states, suggest that a half Fermi gas is realized on the topological surfaces, (b) Spin texture map 
based on spin- ARPES data suggest that the spin-chirality changes sign across the Dirac point, (c) The Dirac node remains well 
defined up a temperature of 300 K suggesting the stability of topological effects up to the room temperature, (d) The Dirac 
cone measured at a temperature of 10 K. (e) Full Dirac cone. [Adapted from D. Hsieh et al, Nature 460, 1101 (2009). [l7(]. 



surface states of Sb, we introduce such a perturbation 
by randomly substituting Bi into the Sb crystal matrix. 
Another motivation for performing such an experiment 
is that the formalism developed by Fu and Kane [lO^ to 
calculate the Z2 topological invariants relies on inversion 
symmetry being present in the bulk crystal, which they 
assumed to hold true even in the random alloy Bii_a^Sba^. 
However, this formalism is simply a device for simplify- 
ing the calculation and the non-trivial = 1 topolog- 
ical class of Bii_a;Sba^ is predicted to hold true even in 
the absence of inversion symmetry in the bulk crystal 
[1^. Therefore introducing light Bi substitutional disor- 
der into the Sb matrix is also a method to examine the ef- 
fects of alloying disorder and possible breakdown of bulk 
inversion symmetry on the surface states of Sb(lll). We 
have performed spin-integrated ARPES measurements 
on single crystals of the random alloy Sbo.gBio.i. Fig- 
ure [13] shows that both the surface band dispersion along 
f -M as well as the surface state Fermi surface retain the 



same form as that observed in Sb(lll), and therefore the 
'topological metal' surface state of Sb(lll) fully survives 
the alloy disorder. Since Bi alloying is seen to only affect 
the band structure of Sb weakly, it is reasonable to as- 
sume that the topological order is preserved between Sb 
and Bio.9iSbo.o9 as we observed. 

In a simpler fashion compared to Bii-^^Sb^^, the topo- 
logical insulator behavior in Bi2Se3 is associated with a 
single band inversion at the surface Brillouin zone center. 
Owing to its larger bandgap compared with Bii_a^Sba^, 
ARPES has shown that its topological properties are pre- 
served at room temperature [l^- Two defining proper- 
ties of topological insulators spin-momentum locking of 
surface states and tt Berry phase can be clearly demon- 
strated in the Bi2Se3 series. The surface states are ex- 
pected to be protected by time-reversal symmetry sym- 
metry, which implies that the surface Dirac node should 
be robust in the presence of nonmagnetic disorder but 
open a gap in the presence of time-reversal symmetry 
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FIG. 15: Protection by time reversal symmetry: Topo- 
logical surface states are robust in the presence of strong non- 
magnetic disorder but open a gap in the presence of time re- 
versal symmetry breaking magnetic impurities and disorder, 
(a) Magnetic impurity such as Fe on the surface of Bi2Se3 
opens a gap at the Dirac point. The magnitude of the gap 
is set by the interaction of Fe ions with the Se surface and 
the time reversal symmetry breaking disorder potential in- 
troduced on the surface, (b) A comparison of surface band 
dispersion with and without Fe doping. (c,d) Non-magnetic 
disorder created via molecular absorbent NO 2 or alkali atom 
adsorption (K or Na) on the surface leaves the Dirac node 
intact in both Bi2Se3 and Bi2Te3. [Adapted from L. Wray et 
aL, Nature Phys. 7, 32 (2011). ^]] 



breaking perturbations. Magnetic impurities such as Fe 
or Mn on the surface of Bi2Se3 open a gap at the Dirac 
point [Figs. [ISfa) and [ISfb)] [22, 'si]. The magnitude 
of the gap is likely set by the interaction of Fe ions with 
the Se surface and the time-reversal symmetry breaking 
disorder potential introduced on the surface. Nonmag- 
netic disorder created via molecular absorbent NO 2 or 
alkali atom adsorption (K or Na) on the surface leaves 



the Dirac node intact 
BisSes and BisTes E 



ts. 15(c) and 15(d)] in both 
These results are consistent 
with the fact that the topological surface states are pro- 
tected by time-reversal symmetry. 

Many of the interesting theoretical proposals that uti- 
lize topological insulator surfaces require the chemical 
potential to lie at or near the surface Dirac point. This 
is similar to the case in graphene, where the chemistry 
of carbon atoms naturally locates the Fermi level at the 
Dirac point. This makes its density of carriers highly 
tunable by an applied electrical field and enables appli- 
cations of graphene to both basic science and microelec- 
tronics. The surface Fermi level of a topological insulator 
depends on the detailed electrostatics of the surface and 
is not necessarily at the Dirac point. Moreover, for nat- 



urally grown Bi2Se3 the bulk Fermi energy is not even 
in the gap. The observed n-type behavior is believed to 
be caused Se vacancies. By appropriate chemical modifi- 
cations, however, the Fermi energy of both the bulk and 
the surface can be controlled. This allowed [17] to reach 
the sweet spot in which the surface Fermi energy is tuned 
to the Dirac point (Fig(T6j). This was achieved by doping 
bulk with a small concentration of Ca, which compen- 
sates the Se vacancies, to place the Fermi level within 
the bulk band gap. The surface was the hole doped by 
exposing the surface to NO2 gas to place the Fermi level 
at the Dirac point. These results collectively show how 
ARPES can be used to study the topological protection 
and tunability properties of the 2D surface of a 3D topo- 
logical insulator. 



FUTURE DIRECTIONS: IDENTIFYING 
MAJORANA PLATFORMS AND TOPOLOGICAL 
SUPERCONDUCTORS 

Recent measurements [25] show that surface instabil- 
ities cause the spin-helical topological insulator band 
structure of Bi2Se3 to remain well defined and non- de- 
generate with bulk electronic states at the Fermi level 
of optimally doped superconducting Cuo.i2Bi2Se3, and 
that this is also likely to be the case for superconduct- 
ing variants of p-type Bi2Te3. These surface states pro- 
vide a highly unusual physical setting in which super- 
conductivity cannot take a conventional form, and is ex- 
pected to realize one of two novel states that have not 
been identified elsewhere in nature. If superconducting 
pairing has even parity, as is nearly universal among the 
known superconducting materials, the surface electrons 
will achieve a 2D non-Abelian superconductor state with 
non-commutative Major ana fermion vortices that can po- 
tentially be manipulated to store quantum information. 
Surface vortices will be found at the end of bulk vortex 
lines as drawn in Fig(T3 If superconducting pairing is 
odd, the resulting state is a novel state of matter known 
as a "topological superconductor" with Bogoliubov sur- 
face quasi-particles present below the superconducting 
critical temperature of 3.8 K. As drawn in FiglTTTc), 
these low temperature surface states would be gapless, 
likely making it impossible to adiabatically manipulate 
surface vortices for quantum computation. The unique 
physics and appHcations of the topological superconduc- 
tor state are distinct from any known material system, 
and will be an exciting vista for theoretical and experi- 
mental exploration if they are achieved for the first time 
in CuajBi2Se3. 
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Tuning the topological helical Dirac electron density 
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FIG. 16: Surface Gating : Tuning the density of helical Dirac electrons to the spin-degenerate Kramers point 
and topological transport regime, (a) A high resolution ARPES mapping of the surface Fermi surface (FS) near F of 
Bi2-(5Ca5Se3 (HI)- The diffuse intensity within the ring originates from the bulk-surface resonance state 11]. (b) The FS after 
0.1 Langmuir (L) of NO2 is dosed, showing that the resonance state is removed, (c) The FS after a 2 L dosage, which achieves 
the Dirac charge neutrality point, (d) High resolution ARPES surface band dispersions through after an NO2 dosage of L, 
(e) 0.01 L, (f) 0.1 L, (g) 0.5 L, (h) 1 L and (i) 2 L. The arrows denote the spin polarization of the bands. We note that due to 
an increasing level of surface disorder with NO 2 adsorption, the measured spectra become progressively more diffuse and the 
total photoemission intensity from the buried Bi2-(5Ca(5Se3 surface is gradually reduced. [Adapted from D. Hsieh et al., Nature 
460, 1101 (2009). 17]]. 
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Growth method for high-quality single crystals 



The undoped and doped Bii_a;Sba; single-crystal sam- 
ples (0 < X < 0.17) used for ARPES experiments were 
each cleaved from a boule grown from a stoichiometric 
mixture of high-purity elements. The boule was cooled 
from 650 °C to 270 °C over a period of five days and was 
annealed for seven days at 270 °C. The samples naturally 
cleaved along the (111) plane, which resulted in shiny 
flat silver surfaces. X-ray diffraction measurements were 
used to check that the samples were single phase, and 
confirmed that the Bio.gSbo.i single crystals presented 
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FIG. 17: A Majorana platform, (a) Topologically protected surface states cross the Fermi level before merging with the 
bulk valence and conduction bands in a lightly doped topological insulator, (b) If the superconducting wavefunction has even 
parity, the surface states will be gapped by the proximity effect, and vortices on the crystal surface will host braidable Majorana 
fermions. (c) If superconducting parity is odd, the material will be a so-called topological superconductor, and new states will 
appear below Tc to span the bulk superconducting gap. (d) Majorana fermion surface vortices are found at the end of bulk 
vortex lines and could be manipulated for quantum computation if superconducting pairing is even. [Adapted from L. Wray 
et al., Nature Phys. 7, 32 (2011). [H] 



in this paper have a rhombohedral A7 crystal structure 
(point group R3m), with room-temperature (T=300K) 
lattice parameters a = 4.51 A and c = 11.78 A indexed us- 
ing a rhombohedral unit cell. The X-ray diffraction pat- 
terns of the cleaved crystals exhibit only the (333), (666), 
and (999) peaks, showing that the cleaved surface is ori- 
ented along the trigonal (111) axis. Room- temperature 
data were recorded on a Bruker D8 diffractometer us- 
ing Cu Ka radiation (A = 1.54 A) and a diffracted-beam 
monochromator. The in-plane crystal orientation was 
determined by Laue X-ray diffraction. During the angle- 
resolved photoemission spectroscopy (ARPES) measure- 
ments a fine alignment was achieved by carefully studying 
the band dispersions and Fermi surface symmetry as an 
internal check for crystal orientation. 



Resistivity characterization 

Temperature-dependent resistivity measurements were 
carried out on single-crystal samples in a Quantum De- 
sign PPMS-9 instrument, using a standard four-probe 
technique on approximately 4x1 xl-mm^, rectangu- 



lar samples with the current in the basal plane, which 
was perpendicular to the trigonal axis. The four con- 
tacts were made by using room-temperature silver paste. 
The data for samples with concentrations ranging from 
x = 0tox = 0.17 showed a systematic change from 
semimetallic to insulating-like behaviour with increasing 
X, in agreement with previously published works p7| , 
which was used as a further check of the antimony con- 
centrations. Conventional magnetic and transport mea- 
surements [l3|, [13, [isl such as these cannot separately 
measure the contributions of the surface and bulk states 
to the total signal. ARPES, on the other hand, is a 
momentum-selective technique [41], which allows for a 
separation of 2D (surface) from 3D (bulk) dispersive en- 
ergy bands. This capability is especially important for 
Bii_a;Sba; because the Dirac point lies at a single point 
in the 3D Brillouin zone, unlike for 2D graphene, where 
the Dirac points can be studied at any arbitrary perpen- 
dicular momentum along a line [39|, [60| . 
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FIG. SI: Method of locating high symmetry bulk reciprocal lattice points of Bio.gSbo.i using incident photon 
energy modulated ARPES. a, Geometry of an ARPES experiment, b, Key parameters relevant to the calculation of the 
positions of the high symmetry points in the 3D BZ. The lattice constants refer to the rhombohedral A7 lattice structure, c, 
Location of L and X points of the bulk BZ in the kx-kz plane together with the constant energy contours that can be accessed 
by changing the angle 0. d, Near Ef intensity map {hu = 55 eV) of the Fermi surface formed by the surface states covering 
an entire surface BZ, used to help locate various in-plane momenta, in units of the photoelectron emission angle along two 
orthogonal spatial directions. The electron pockets near M in FiglSj: (main text) appear as lines in FigjSld due to relaxed 
k-resolution in order to cover a large k-space in a single shot. [Adapted from D. Hsieh et al., Nature 452, 970 (2008) [14]]. 



ARPES 

Spin-integrated angle-resolved photoemission spec- 
troscopy (ARPES) measurements were performed with 
14 to 30 eV photons on beam line 5-4 at the Stanford 
Synchrotron Radiation Laboratory, and with 28 to 32 eV 
photons on beam line 12 at the Advanced Light Source, 
both endstations being equipped with a Scienta hemi- 
spherical electron analyzer (see VG Scienta manufacturer 
website for instrument specifications). Spin-resolved 
ARPES measurements were performed at the SIS beam 
line at the Swiss Light Source using the COPHEE spec- 
trometer ([56] p. 15) with a single 40 kV classical Mott 
detector and photon energies of 20 and 22 eV. The typ- 
ical energy and momentum resolution was 15 meV and 
1.5% of the surface Brillouin zone (BZ) respectively at 
beam line 5-4, 9 meV and 1% of the surface BZ respec- 
tively at beam line 12, and 80 meV and 3% of the sur- 
face BZ respectively at SIS using a pass energy of 3 eV. 
Cleaving these samples in situ between 10 K and 55 K 
at chamber pressures less than 5 xlO~^^ torr resulted in 
shiny flat surfaces, characterized in situ by low energy 



electron diffraction (LEED) to be clean and well ordered 
with the same symmetry as the bulk. This is consistent 
with photoelectron diffraction measurements that show 
no substantial structural relaxation of the Sb(lll) sur- 
face [57]. 



Systematic methods for separating bulk from 
surface electronic states 



ARPES is a photon-in, electron-out technique |41|. 
Photoelectrons ejected from a sample by a monochro- 
matic beam of radiation are collected by a detector capa- 
ble of measuring its kinetic energy Ekin- By varying the 
detector angles, 6 and relative to the sample surface 
normal, the momentum of the photoelectrons, K, can 
also be determined (as illustrated in ISip. By employing 
the commonly used free-electron final state approxima- 
tion, we can fully convert from the measured kinetic en- 
ergy and momentum values of the photoelectron to the 
binding energy, Eb^ and Bloch momentum values k of 
its initial state inside the crystal, via 
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Tight binding bandstructure near L point 
w/ spin-orbit coupling 
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FIG. S2: Spin-orbit coupling has a profound effect on the band structure of bismuth near L point, a, Schematic 
of the bulk 3D BZ and the projected BZ of the (111) surface. b,c, Calculated tight binding band structure of bismuth including 
a spin-orbit coupling strength of 1.5 eV along two orthogonal cuts through the L point in the 1** bulk BZ. d,e, Tight binding 
band structure along the same two directions as b and c calculated without spin-orbit coupling. The inter-band gap of 13.7 
meV is barely visible on the scale of b and c. [Adapted from D. Hsieh et aL, Nature 452, 970 (2008) [13]]. 

from 17 eV to 55 eV, energy resolutions ranging from 9 
meV to 40 meV and momentum resolution better than 
1.5% of the surface Brillouin zone, using Scienta elec- 
tron analysers. The combination of high spatial resolu- 
tion and high crystalline quality enabled us to probe only 
the highly ordered and cleanest regions of our samples. 
Single-crystal Bii_a;Sba^ samples were cleaved in situ at 
a temperature of 15 K and chamber pressures less than 
8 X 10~^^ torr, and high surface quality was checked 
throughout the measurement process by monitoring the 
EDC linewidths of the surface state. To measure the 
near-£^i? dispersion of an electronic band along a direc- 
tion normal to the sample surface, such as the direction 
from X(27r/V^a,0,87r/c) to L(27r/V3a, 0, IItt/c) shown 
in Fig. [3^, EDCs were taken at several incident photon 
energies. The kinetic energy of the photoelectron at Ef 
is different for each value of hu, so the angle was first 
adjusted and then held fixed for each hu so as to keep 
kx constant at 2ixl\fZa = 0.8 A~^ for electrons emitted 
near Ep- To ensure that the in-plane momentum re- 
mained constant at M (the L-X line projects onto M) 



\EB\ = hv-W -E^,n 
kx = = ^ \/ 2meEkin^\nO 

kz = ^\/2me{EkinCos'^0 - Vo) 

where we have set cp = 0^ W is the work function, mg is 
the electron mass and Vq is an experimentally determined 
parameter, which is approximately — 10 eV for bismuth 
|6ll , l62| . Features in the ARPES spectra originating from 
bulk initial states (dispersive along the /c^-direction) were 
distinguished from those originating from surface initial 
states (non-dispersive along the /c^-direction) by study- 
ing their dependence on incident photon energy, hu, and 
converting this to dependence on kz via the displayed 
equations. ARPES data were collected at beamlines 
12.0.1 and 10.0.1 of the Advanced Light Source at the 
Lawrence Berkeley National Laboratory, as well as at 
the PGM beamline of the Synchrotron Radiation Cen- 
ter in Wisconsin, with incident photon energies ranging 
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for each EDC, a complete neaT-Ep intensity map was 
generated for each photon energy to precisely locate the 
M-point (see FiglSTH). We note that because the bulk 
crystal has only three-fold rotational symmetry about the 
/c^-axis, the reciprocal lattice does not have mirror sym- 
metry about the kx = plane. Therefore, scans taken 
at -\-6 and -6 for the same photon energy probe different 
points in the bulk 3D Brillouin zone; this is responsible 
for the absence of the bulk A-shaped band in Fig. |4f. 



Confirming the bulk nature of electronic bands by 
comparison with theoretical calculations 

In an ARPES experiment fFiglSlb). three dimensional 
(3D) dispersive bulk electronic states can be identified as 
those that disperse with incident photon energy, whereas 
surface states do not. As an additional check that 
we have indeed correctly identified the bulk bands of 
Bio.gSbo.i in Figs [2] and [3l we also measured the disper- 
sion of the deeper lying bands well below the Fermi level 
(Ep) and compared them to tight binding theoretical 
calculations of the bulk bands of pure bismuth following 
the model of Liu and Allen (1995) [25]. A tight-binding 
approach is known to be valid since Bio.gSbo.i is not a 
strongly correlated electron system. As Bio.gSbo.i is a 
random alloy (Sb does not form a superlattice [24]) with 
a relatively small Sb concentration (~0.2 Sb atoms per 
rhombohedral unit cell), the deeper lying band structure 
of Bio.gSbo.i is expected to follow that of pure Bi be- 
cause the deeper lying (localized wave function) bands 
of Bio.gSbo.i are not greatly affected by the substitu- 
tional disorder, and no additional back folded bands are 
expected to arise. Since these deeper lying bands are 
predicted to change dramatically with /c^, they help us 
to finely determine the experimentally probed kz values. 
Fig|S2lF shows the ARPES second derivative image (SDI) 
of a cut parallel to KMK that passes through the L point 
of the 3D Brillouin zone (BZ), and FiglS2h shows a par- 
allel cut that passes through the 0.3 XL point (FiglS2b). 
The locations of these two cuts in the 3D bulk BZ were 
calculated from the ARPES kinematic relations, from 
which we can construct the constant energy contours 
shown in Fig. [STb . By adjusting such that the in- 
plane momentum kx is fixed at approximately 0.8 A~^ 
(the surface M point), at a photon energy hu =29 eV, 
electrons at the Fermi energy {Eb=0 eV) have a kz that 
corresponds to the L point in the 3^^ bulk BZ. By ad- 
justing such that the in-plane momentum kx is fixed at 
approximately -0.8 A~^, at a photon energy hu = 20 eV, 
electrons at a binding energy of -2 eV have a kz near 0.3 
XL. 

There is a clear kz dependence of the dispersion of mea- 
sured bands A, B and C, pointing to their bulk nature. 
The bulk origin of bands A, B and C is confirmed by their 
good agreement with tight binding calculations (bands 3, 
4 and 5 in Figs [S2^ and i) , which include a strong spin- 
orbit coupling constant of 1.5 eV derived from bismuth 



(25j . Band 3 drops below -5 eV at the 0.3 XL point. 
The slight differences between the experimentally mea- 
sured band energies and the calculated band energies at 
ky = A~^ shown in Figl52t-i are due to the fact that 
the ARPES data were collected in a single shot, taken 
in constant mode. This means that electrons detected 
at different binding energies will have slightly different 
values of kz, whereas the presented tight binding calcu- 
lations show all bands at a single kz. We checked that 
the magnitude of these band energy differences is indeed 
accounted for by this explanation. Even though the La 
and Lg bands in Bio.gSbo.i are inverted relative to those 
of pure semimetallic Bi, calculations show that near Ep^ 
apart from an insulating gap, they are "mirror" bands in 
terms of k dispersion (see bands 5 and 6 in FigIS2^). Such 
a close match to calculations, which also predict a linear 
dispersion along the kz cut near Ep f FiglSTH). provides 
strong support that the dispersion of band C, near Ep^ is 
in fact linear along kz. Focusing on the A-shaped valence 
band at L, the EDCs fFiglS2k) show a single peak out to 
ky ~ ±0.15 A~^ demonstrating that it is composed of a 
single band feature. Outside this range however, an addi- 
tional feature develops on the low binding energy side of 
the main peak in the EDCs, which shows up as two well 
separated bands in the SDI image (Figl3f) and signals a 
splitting of the band into bulk representative and surface 
representative components (Figl52k.f). Unlike the main 
peak that disperses strongly with incident photon energy, 
this shoulder-like feature is present and retains the same 
A-shaped dispersion near this k-region (open circles in 
Figs [S2^ and i) for all photon energies used, supporting 
its 2D surface character. This behaviour is quite unlike 
bulk band C, which attains the A-shaped dispersion only 
near 29 eV (see main text Fig. (SJd). 



Spin-orbit coupling is responsible for the unique 
Dirac-like dispersion behaviour of the bulk bands 
near Ef 

According to theoretical models, a strongly spin-orbit 
coupled bulk band structure is necessary for topologi- 
cal surface states to exist Therefore it 
is important to show that our experimentally measured 
bulk band structure of Bio.gSbo.i can only be accounted 
for by calculations that explicitly include a large spin- 
orbit coupling term. As shown in the previous section, 
the measured bulk band dispersion of Bio.gSbo.i gener- 
ally follows the calculated bulk bands of pure Bi from 
a tight binding model. The dispersion of the bulk va- 
lence and conduction bands of pure bismuth near Ep dit 
the L point from such a tight binding calculation [25| 
with a spin-orbit coupling constant of 1.5 eV are shown 
in Fig. [S3b and c, which show a high degree of linear- 
ity. The high degree of linearity can be understood from 
a combination of the large Fermi velocity (vi? ~ 6 eV 
A along ky) and small inter-band (below Ep) gap A = 
13.7 meV (Fig. [S3|) . This calculated inter-band gap of 
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FIG. S3: Identification of the bulk band features of Bio.gSbo.i. Experimental band-structure determined by 
ARPES is compared to bulk tight-binding calculations of bismuth to further identify the deeper lying bulk 
bands and their symmetry origins, a, Energy distribution curves (EDCs) along a k-space cut given by the upper yellow 
line shown in schematic c which goes through the bulk L point in the 3^^ BZ (hu — 29 eV). The corresponding ARPES intensity 
in the vicinity of L is shown in e. b, EDCs along the lower yellow line of c which goes through the point a fraction 0.3 of the 
k-distance from A to L (hu = 20 eV), showing a dramatic change of the deeper lying band dispersions. (This cut was taken 
at a kx value equal in magnitude but opposite in sign to that in a as described in the SM text). f,h, The ARPES second 
derivative images (SDI) of the raw data shown in a and b to reveal the band dispersions. The flat band of intensity at Ef is 
an artifact of taking SDI. g,i, Tight binding band calculations of bismuth including spin-orbit coupling, using Liu and Allen 
model [25], along the corresponding experimental cut directions shown in f and h. The bands (colored solid lines) labelled 3 to 
6 are derived from the symmetries associated with the 6p-orbitals and their dispersion is thus strongly influenced by spin-orbit 
coupling. The inter-band gap between bands 5 and 6 is barely visible on the scale of Fig. [S2k . The circled curves mark the 
surface state dispersion, which is present at all measured photon energies (no kz dispersion). There is a close match of the bulk 
band dispersion between the data and calculations, confirming the presence of strong spin-orbit coupling, d, Tight binding 
valence band (5) dispersion of bismuth in the ky-kz momentum plane showing linearity along both directions. The close match 
between data and calculation along ky suggests that the dispersion near Ef along kz is also linear. [Adapted from D. Hsieh et 
al., Nature 452, 970 (2008) [14]]. 



Bi (13.7 meV) is smaller than our measured lower limit 
of 50 meV (main text Fig. [2^) for the insulating gap 
of Bio.gSbo.i. To illustrate the importance of spin-orbit 
coupling in determining the band structure near L, we 
show the dispersion along ky and kz calculated without 
spin-orbit coupling (Fig. \S3\ 1 and e). While the disper- 
sion along kz is not drastically altered by neglecting the 
spin-orbit coupling, the dispersion along ky changes from 
being linear to highly parabolic. This is further evidence 
that our measured Dirac point can be accounted for only 
by including spin-orbit coupling. A strong spin- orbit cou- 



pling constant acts as an internal quantizing magnetic 
field for the electron system [32] which can give rise to a 
quantum spin Hall effect without any externally applied 
magnetic field ,65^ ? , i64]- Therefore, the existence 
or the spontaneous emergence of the surface or boundary 
states does not require an external magnetic field. 
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FIG. S4: The Kramers' point, the gapless nature and topology of surface states in insulating Bio.gSbo.i is revealed 
through high spatial and k-resolution ARPES. a, Energy distribution curves (EDCs) of a low resolution ARPES scan 
along the surface F-M cut of Bio.gSbo.i. b, The surface band dispersion second derivative image along F-M obtained by piecing 
together four high resolution ARPES scans. See main text Fig|4]for explanation of other features, c, EDCs of high resolution 
ARPES scans in the vicinity of surface Fermi crossings 1 and 2 and crossings 3, 4 and 5 (left panels). These crossings form the 
surface Fermi surface shown in the upper right panel of c (see also main text Fig|3]). High resolution ARPES scans along cut 
directions A, B and C are further evidence for a surface Fermi surface. [Adapted from D. Hsieh et al., Nature 452, 970 (2008) 

EJ 



Matching the surface state Fermi crossings and the 
topology of the surface Fermi surface in bulk 
insulating Bio.gSbo.i 

In order to count the number of singly degenerate sur- 
face state Fermi crossings [l, |42|, [g^ along the f -M cut of 
the surface BZ, high photon energy ARPES scans, which 
allow mapping of the entire k range from f-M to fall 
within the detector window at the expense of lower in- 
strument resolution, were taken to preliminarily identify 
the k-space locations of the Fermi crossings (Fig. [S4k). 
Having determined where these surface state Fermi cross- 
ings lie in k-space, we performed several high resolution 
ARPES scans, each covering a successive small k inter- 
val in the detector window, in order to construct a high 
resolution band mapping of the surface states from f to 
M. The second derivative image of the surface band dis- 
persion shown in FiglS4b was constructed by piecing to- 



gether four such high resolution scans. FiglS4h shows en- 
ergy distribution curves of high resolution ARPES scans 
in the vicinity of each surface Fermi crossing, which to- 
gether give rise to the surface Fermi surface shown. No 
previous work jsl, Isl - lssf has reported the band disper- 
sion near the L-point (thus missing the Dirac bands) or 
resolved the Kramers point near the M point, which is 
crucial to determine the topology of the surface states. 
For this reason there is no basis for one-to-one compari- 
son with previous work, since no previous ARPES data 
exists in the analogous k-range. Note that surface band 
dispersions along the cuts A, B and C are highly lin- 
ear. This is indirect evidence for the existence of the 
bulk Dirac point since surface states are formed when the 
bulk state wave functions are subjected to the boundary 
conditions at the cleaved plane. 
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Two-step fitting analysis procedure of Spin-Resolved 
ARPES measurements of insulating Bii-^Sbx 

Here we present details of the spin-resolved ARPES 
analysis on bulk insulating Bio.91Sbo.09 that show how 
we arrive at a spin-resolved surface band dispersion such 
as that presented in Figure [5jG) in the main text. In the 
VUV incident photon energy regime that we use, spin 
conserving photoemission processes (where the electric 
field of light only acts on the orbital degree of freedom 
of the electron inside a solid) dominate over spin non- 
conserving processes (which arise from coupling to the 
magnetic field of light) [ss']. Therefore we are confident 
that the photo-emission process does not change the spin 
polarization of the electrons. Figure [S5f B) shows a spin 
averaged momentum distribution curve (MDC) along the 
r to -M direction taken dit Eb = -25 meV, indicated by 
the green line shown in Figure ISSf A). This MDC was 
obtained by summing the signal coming from both left 
and right electron detectors in the Mott polarimeter (see 
diagram in Fig 19]^ A) of the main text). Lorentzian line- 
shapes denoted P and a non-polarized background B are 
fitted to this MDC, which are used as inputs to the two- 
step fitting routine developed by Meier et al. [45] in the 
following way. To begin with, a spin polarization vector 
= {PI' ^ ^1' ^ PI') = (cos 6>i COS ^i, COS 6>i sin ^i, sin l9i) is 
assigned to each band, where Oi and are referenced 
to the primed Mott coordinate frame. Here it is neces- 
sary to assume a spin magnitude of one because only two 
spin components are measured by a single Mott detec- 
tor. Such an assumption is likely valid since even though 
the spin polarization is no longer a good quantum num- 
ber due to spin-orbit coupling, the bands near f are ex- 
pected to exhibit a high degree of spin polarization since 
the spin-orbit coupling is smallest near f . Moreover, 
common strong spin-orbit coupled materials such as gold 
have been experimentally shown to exhibit 100% spin 
polarized surface states [46l] . A spin-resolved spectrum is 
then defined for each peak i using /^'^'^ = /*(! ± P^)/6, 
where a = x'^y'^z'^ and + and — correspond to the 
spin direction being parallel (t) or antiparallel (|) to 
a. The full spin-resolved spectrum is then given by 
/^'^ = ^a^'^ + where B is the unpolarized back- 
ground, from which the spin polarization of each spatial 
component can be obtained as = (/^ — /^)/(/^ + /^). 
This latter expression is a function of Oi and and is 
used to fit to the experimental data. 

The spin polarization data for the y' and z' compo- 
nents (i.e. Py' and P^/) are obtained by taking the dif- 
ference between the intensities of the left-right (or top- 
bottom) electron detectors over their sum, all divided 
by the Sherman function, which is calibrated using the 
methods in ([56], p. 36). Typical electron counts on the 
detector reach 5 x 10^, which places an error bar of ap- 
proximately ±0.01 for each point on our polarization 
curves. To account for unequal sensitivities between a 
detector pair, we applied a small multiplicative factor to 
the intensity from one detector to ensure that the unpo- 



larized background intensity yields zero polarization. Re- 
sultant curves are shown in Figure [S5jC). The best fit pa- 
rameters (P^/, P^/, P^/), which are expressed in the sam- 
ple coordinates through an appropriate coordinate trans- 
formation [inset of Fig. lSSl C)] are shown in Figure lSSl D). 
Even though the measured polarization only reaches a 
magnitude of around ±0.1, this is similarly seen in stud- 
ies of Bi thin films [42| and is due to the non-polarized 
background and overlap of adjacent peaks with different 
spin polarization. These effects are extremely sensitive 
to the sample alignment due to the very narrow Fermi 
surface features. The fitted parameters [Fig. [S5l D)] are 
consistent with spins being nearly aligned along the ±^ 
direction, with bands 11 and rl having nearly opposite 
spin as required by time reversal symmetry, and with 
these spins nearly parallel to those of 11 and rl respec- 
tively measured for Sb [main text FiglHl^F)]. The small 
departures from ideality likely originate from the scan 
direction not being exactly along f-M. Bands 11 and 12 
display opposite spin, which indicates that they form a 
Kramers pair split by spin-orbit coupling, and the fact 
that bands 12 and /3 have the same spin suggests that 
they originate from the same band. 

To show that bands 12 and /3 connect above Ep as we 
have drawn in Figure [S5f A), and are thus expected to 
have the same spin, we map the surface band dispersion 
of Te doped Bii_a;Sba^ that is known to be an electron 
donor [66|. Figure lSSf F) shows that the hole band formed 
by crossings 2 and 3 in insulating Bii-a^Sb^^ [Fig. [55f A)] 
has sunk completely below Ep with 0.5% Te doping, and 
is in fact the same band. 



Method of using incident photon energy modulated 
ARPES to separate the bulk from surface electronic 
states of Sb 

In this section we detail incident photon energy mod- 
ulated ARPES experiments on the low lying electronic 
states of single crystal Sb(lll), which we employ to iso- 
late the surface from bulk-like electronic bands over the 
entire BZ. Figure [S6F C) shows momentum distributions 
curves (MDCs) of electrons emitted at Ep as a function 
oi kx (II f-M) for Sb(lll). The out-of-plane component 
of the momentum kz was calculated for different inci- 
dent photon energies {hv) using the free electron final 
state approximation with an experimentally determined 
inner potential of 14.5 eV (HH, [69|. There are four peaks 
in the MDCs centered about V that show no dispersion 
along kz and have narrow widths of /S.kx ^ 0.03 A~^. 
These are attributed to surface states and are similar to 
those that appear in Sb(lll) thin films [68]. As hv is 
increased beyond 20 eV, a broad peak appears at kx ^ 
-0.2 A~^, outside the k range of the surface states near 
r, and eventually splits into two peaks. Such a strong kz 
dispersion, together with a broadened linewidth (A/c^^ ^ 
0.12 A~^), is indicative of bulk band behavior, and indeed 
these MDC peaks trace out a Fermi surface [Fig. [S6fD)] 
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FIG. S5: (A) The surface band dispersion ARPES second derivative image (SDI) along the T to -M direction of bulk insulating 
Bio.9iSbo.o9- Dashed white lines are guides to the eye. The intensity of bands /4, 5 is scaled up for clarity. (B) MDC of the 
spin averaged spectrum at Eb = -25 meV [green line below Ef in (A)] using a photon energy hu = 22 eV, together with the 
Lorentzian peaks of the fit. (C) Measured spin polarization curves (symbols) for the y' and z' (Mott coordinate) components 
together with the fitted lines (SM). The relative orientation of the sample (un-primed) to Mott (primed) coordinates is shown 
in the inset. The polar angle is rotated during the measurement to access different values of kx. At normal emission (^ = 0°), 
the z' and z axes are parallel and the y' axis is rotated from the y axis by 45°. (D) The in-plane and out-of-plane spin 
polarization components in the sample coordinate frame obtained from the spin polarization fit. The symbols refer to those 
in (B). The fitted parameters are consistent with 100% polarized spins. (E) Spin resolved spectra for the y component based 
on the fitted spin polarization curves shown in (C). Spin up (down) refers to the spin direction being approximately parallel 
to the +(-)y direction. (F) The surface band dispersion SDI centered about F of (Bio. 925 Sbo. 075)0. 995Teo.oo5- Electron doping 
through Te reveals that bands 12 and /3 are connected above Ep. [Adapted from D. Hsieh et a/., Science 323, 919 (2009) [3]. 



that is similar in shape to the hole pocket calculated for 
bulk Sb near H [67^ . Therefore by choosing an appropri- 
ate photon energy (e.g. < 20 eV), the ARPES spectrum 
8it Ef along F-M will have contributions from only the 
surface states. The small bulk electron pocket centered 
at L is not accessed using the photon energy range we 
employed [Fig. [HUD)]. 

Now we describe the experimental procedure used to 
distinguish pure surface states from resonant states on 
Sb(lll) through their spectral signatures. ARPES spec- 
tra along f -M taken at three different photon energies are 
shown in Fig. [S7l Near f there are two rather linearly 
dispersive electron like bands that meet exactly at f at a 
binding energy Eb ~ -0.2 eV. This behavior is consistent 
with a pair of spin-split surface bands that become de- 
generate at the time reversal invariant momentum (A:t) 
F due to Kramers degeneracy. The surface origin of this 
pair of bands is established by their lack of dependence 
on hv [Fig. [S7lf A)-(C)]. A strongly photon energy disper- 
sive hole like band is clearly seen on the negative kx side 
of the surface Kramers pair, which crosses Ef for hu = 
24 eV and gives rise to the bulk hole Fermi surface near 
H [Fig. [SgD)]. For hu < 20 eV, this band shows clear 



back folding near Eb ^ -0.2 eV indicating that it has 
completely sunk below Ef- Further evidence for its bulk 
origin comes from its close match to band calculations 
[Fig. IS6lfD)]. Interestingly, at photon energies such as 18 
eV where the bulk bands are far below Ef^ there re- 
mains a uniform envelope of weak spectral intensity near 
Ef in the shape of the bulk hole pocket seen with hu = 
24 eV photons, which is symmetric about f . This enve- 
lope does not change shape with hu suggesting that it 
is of surface origin. Due to its weak intensity relative to 
states at higher binding energy, these features cannot be 
easily seen in the energy distribution curves (EDCs) in 
Fig. A)- (C), but can be clearly observed in the MDCs 
shown in Fig. [S6lf C) especially on the positive kx side. 
Centered about the M point, we also observe a crescent 
shaped envelope of weak intensity that does not disperse 
with kz [Fig. [S7F D)-(F)]. pointing to its surface origin. 
Unlike the sharp surface states near f , the peaks in the 
EDCs of the feature near M are much broader {AE ~80 
meV) than the spectrometer resolution (15 meV). The 
origin of this diffuse ARPES signal is not due to sur- 
face structural disorder because if that were the case, 
electrons at f should be even more severely scattered 
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FIG. S6: (A) Schematic of the bulk BZ of Sb and its (111) 
surface BZ. The shaded region denotes the momentum plane 
in which the following ARPES spectra were measured. (B) 
LEED image of the in situ cleaved (111) surface exhibiting 
a hexagonal symmetry. (C) Select MDCs at taken with 
photon energies from 14 eV to 26 eV in steps of 2 eV, taken 
in the TXLU momentum plane. Peak positions in the MDCs 
were determined by fitting to Lorentzians (green curves). (D) 
Experimental 3D bulk Fermi surface near H (red circles) and 
2D surface Fermi surface near F (open circles) projected onto 
the kx-kz plane, constructed from the peak positions found in 
(C). The kz values are determined using calculated constant 
hiy contours (black curves) (SM). The shaded gray region is 
the theoretical hole Fermi surface calculated in |67]. [Adapted 
from D. Hsieh et a/., Science 323, 919 (2009) Hi]]. 



from defects than those at M. In fact, the occurrence of 
both sharp and diffuse surface states originates from a 
k dependent coupling to the bulk. As seen in Fig|7jD) 
of the main text, the spin-split Kramers pair near f lie 
completely within the gap of the projected bulk bands 
near Ef attesting to their purely surface character. In 
contrast, the weak diffuse hole like band centered near 
kx = 0.3 A~^ and electron like band centered near kx 
= 0.8 A~^ lie completely within the projected bulk va- 
lence and conduction bands respectively, and thus their 
ARPES spectra exhibit the expected lifetime broadening 
due to coupling with the underlying bulk continuum ^] . 



Method of counting spin Fermi surface kr enclosures 
in pure Sb 

In this section we give a detailed explanation of why 
the surface Fermi contours of Sb(lll) that overlap with 
the projected bulk Fermi surfaces can be neglected when 
determining the uq class of the material. Although the 
Fermi surface formed by the surface resonance near M en- 
closes the /ct M, we will show that this Fermi surface will 
only contribute an even number of enclosures and thus 
not alter the overall evenness or oddness of kr enclosures. 
Consider some time reversal symmetric perturbation that 
lifts the bulk conduction band completely above Ef 
so that there is a direct excitation gap at L. Since this 
perturbation preserves the energy ordering of the and 
Ls states, it does not change the class. At the same 
time, the weakly surface bound electrons at M can evolve 
in one of two ways. In one case, this surface band can also 
be pushed up in energy by the perturbation such that it 
remains completely inside the projected bulk conduction 
band [Fig. ISSlfA)]. In this case there is no more density 
of states dit Ef around M. Alternatively the surface band 
can remain below Ef so as to form a pure surface state 
residing in the projected bulk gap. However by Kramers 
theorem, this SS must be doubly spin degenerate at M 
and its FS must therefore enclose M twice [Fig. lSSf B)]. In 
determining for semi- metallic Sb(lll), one can there- 
fore neglect all segments of the FS that lie within the 
projected areas of the bulk FS [Fig|7fG) of main text] 
because they can only contribute an even number of FS 
enclosures, which does not change the modulo 2 sum of 
kr enclosures. 

In order to further experimentally confirm the topolog- 
ically non-trivial surface band dispersion shown in fig- 
ures EJC) and (D) of the main text, we show ARPES 
intensity maps of Sb(lll) along the -K— f — K direction. 
Figure ISSlfC) shows that the inner V-shaped band that 
was observed along the -M— f — M direction retains its 
V-shape along the -K— f — K direction and continues to 
cross the Fermi level, which is expected since it forms 
the central hexagonal Fermi surface. On the other hand, 
the outer V-shaped band that was observed along the - 
M— f — M direction no longer crosses the Fermi level along 
the -K— f — K direction, instead folding back below the 
Fermi level around ky = 0.1 A~^ and merging with the 
bulk valence band [Fig. IBSF C)]. This confirms that it is 
the Si (2) band starting from F that connects to the bulk 
valence (conduction) band, in agreement with the calcu- 
lations shown in figure [Tl^D) of the main text. 

Physical interpretation of um : the mirror Chern 
number and an analogy with the spin-Chern number 

In this section we will describe how a mirror Chern 
number arises from the crystal symmetry of Bii-a^Sb^^. 
Electronic states in the mirror plane {ky = 0) [Fig.lSOlfA)] 
are eigenstates of the mirror operator M{y) with eigen- 
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FIG. S7: ARPES intensity maps of Sb(lll) as a function of kx near f (A)-(C) and M (D)-(F) and their corresponding energy 
distribution curves, taken using hu = 24 eV, 20 eV and 18 eV photons. The intensity scale of (D)-(F) is a factor of about 
twenty smaller than that of (A)-(C) due to the intrinsic weakness of the ARPES signal near M. [Adapted from D. Hsieh et 
aL, Science 323, 919 (2009) 0]. 



values ±i. M{y) is closely related to, but not exactly 
the same as the spin operator Sy. It may be written 
as M{y) = PC2{y)' the product of the parity operator 
P : {x^y^z) (—X, —y^ —z) and a twofold rotation op- 
erator C2{y)' (x^y^z) {—x^y^—z). For a free spin, 
P does not affect the pseudovector spin, and C2{y) sim- 
ply rotates the spin. Thus, M{y) = ex.-p[—i7TSy/h]. For 
spin eigenstates Sy = ±h/2^ this gives M{y) = ^i. In 
the crystal with spin-orbit interaction on the other hand, 
Sy is no longer a good quantum number, but M{y) still 
is. The energy bands near the Fermi energy in Bii-a^Sba^ 
are derived from states with even orbital mirror symme- 
try and satisfy M{y) oc —i sign{{Sy))^ as detailed in (ioj 
and summarized below. 

Unlike the bulk states which are doubly spin degen- 
erate, the surface state spin degeneracy is lifted due to 
the loss of crystal inversion symmetry at the surface, giv- 
ing rise to the typical Dirac like dispersion relations near 
time reversal invariant momenta [Fig. [59F B)fc(C)]. For 
surface states in the mirror plane ky = with M{y) = ±z, 
the spin split dispersion near kx = has the form 
E = zthvkx. Assuming no other band crossings occur, 
the sign of the velocity v is determined by the topologi- 
cal mirror Chern number {um) describing the bulk band 
structure. When um = 1, the situation in figure [59f B) 



is realized where it is the spin up {{Sy) \\ y) band that 
connects the bulk valence to conduction band going in 
the positive kx direction (i.e. the spin up band has a 
velocity in the positive x direction). For um = — 1 the 
opposite holds true [Fig. ISOF C)]. These two possibilities 
also lead to two distinct chiralities of the central Fermi 
surface as shown in figures [S9lfB)fc(C). From our spin- 
resolved ARPES data on both insulating Bii-a^Sb^^ and 
pure Sb, we find that the surface polarized band disper- 
sions are consistent with um = — 1 [Figs [S9F D)fc(E)]. 
suggesting that their bulk electron wavefunctions exhibit 
the anomalous value um = — 1 predicted in (42), which is 
not realizable in free electron systems with full rotational 
symmetry. 

There is an intimate physical connection between a 
2D quantum spin Hall insulator and the 2D k-space mir- 
ror plane of a 3D topological insulator. In the former 
case, the occupied energy bands for each spin eigenvalue 
will be associated with an ordinary Chern integer 
from which a non-zero spin- Chern number can be defined 
Us = {n^—ni)/2. In the latter case, it is the mirror eigen- 
value of the occupied energy bands that have associated 
with them Chern integers n+i^_i, from which a non-zero 
mirror Chern number can be defined um = (n+i— n_i)/2. 
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FIG. S8: (A) Schematic of the surface band structure of 
Sb(lll) under a time reversal symmetric perturbation that 
hfts the bulk conduction (La) band above the Fermi level 
(Ef)- Here the surface bands near M are also lifted com- 
pleted above Ef- (B) Alternatively the surface band near M 
can remain below Ef itl which case it must be doubly spin 
degenerate at M. (C) ARPES intensity plot of the surface 
states along the -K— F— K direction. The shaded green regions 
denote the theoretical projection of the bulk valence bands, 
calculated using the full potential linearized augmented plane 
wave method using the local density approximation including 
the spin-orbit interaction (method described in 40). Along 
this direction, it is clear that the outer V-shaped surface band 
that was observed along the -M— F— M now merges with the 
bulk valence band. [Adapted from D. Hsieh et al., Science 
323, 919 (2009) il6[]. 
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